Casimir forces between arbitrary compact objects 
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We develop an exact method for computing the Casimir energy between arbitrary compact objects, 
either dielectrics or perfect conductors. The energy is obtained as an interaction between multipoles, 
generated by quantum current fluctuations. The objects' shape and composition enter only through 
their scattering matrices. The result is exact when all multipoles are included, and converges 
rapidly. A low frequency expansion yields the energy as a series in the ratio of the objects' size 
to their separation. As an example, we obtain this series for two dielectric spheres and the full 
interaction at all separations for perfectly conducting spheres. 
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The electromagnetic (EM) force between neutral bod- 
ies is governed by the coordinated dance of fluctuating 
charges pQ. At the atomic scale, this attractive interac- 
tion appears in the guises of van der Waals, Keesom, De- 
bye, and London forces. The collective behavior of con- 
densed atoms is better formulated in terms of dielectric 
properties. In 1948, Casimir computed the force between 
two parallel metallic plates by focusing on the quantum 
fluctuations of the EM field between the two plates [2]. 
This was extended by Lifshitz to dielectric plates, ac- 
counting for the fluctuating fields in the media [3 . The 
force between atoms at asymptotically large distances 
was computed by Casimir and Polder [1] and related 
to the atoms' polarizabilities. For compact objects, such 
as two spheres, Feinberg and Sucher [5] generalized this 
work to include magnetic effects. 

In this Letter we obtain the EM Casimir interaction 
between compact objects at arbitrary separations [BJ, and 
determine explicitly the dependence on shape and mate- 
rial properties [7]. In a qualitative sense, our approach 
is similar to a multipole expansion for the fluctuating 
sources. The dependence on shape and material appears 
through the susceptibility to current fluctuations, and 
is related to the scattering of EM waves by the object. 
While the scattering matrix is in principle complicated, 
there are tools for computing it and it is known for cer- 
tain geometries. As an example, we compute the EM 
force between two dielectric spheres at any separation. 

Earlier studies of the Casimir force between compact 
objects include a multiple reflection formalism[8j, which 
in principle could be applied to perfect conductors of ar- 
bitrary shape. A formulation of the Casimir energy of 
compact objects in terms of their scattering matrices, for 
a scalar field coupled to a dielectric background, is intro- 
duced in Ref. [9], where it is suggested that it can also 
be extended to the EM case. 



Many of our results can be derived by either Green's 
function or path integral methods. We shall sketch the 
latter derivation — due to the letter format only the key 
steps are outlined, and details are left for a more com- 
plete exposition [10] . Note first that since the objects are 
fixed in time, the action is diagonal in the frequency k. 
Therefore in all subsequent steps we can treat each fre- 
quency independently, and integrate over k at the end. 
The Casimir energy can be associated with modifications 
of gauge field fluctuations due to constraints imposed by 
boundary conditions at the material objects. An alter- 
native and equivalent description, stressed by Schwingcr 
, is to attribute the Casimir interaction to fluctuating 
current and charge densities J, g inside the objects. In 
the latter formulation, the EM gauge and scalar potential 
[A(x, t),<&(x, t)] are given for each source configuration 
by the classical solutions, which in Lorcntz gauge read 

[A(x),$(x)] = J dx'G (x,x')[J(x'),e(x')] , (1) 

with Go(x,x') = e lfc l x_x 'l/(47r|x — x'|). For path integral 
quantization, we integrate over all allowed configurations 
of the fluctuating currents, weighted by the appropriate 
action. The Lagrangian for a collection of currents in 
vacuum is the kinetic energy ^JA minus the potential 
energy \g<&- This yields, using Eq. and the continuity 
equation VJ = ikg, the action S[3] = J(dfc/47r)(Sfc[J] + 
<Sfc[J]) for the current densities {3 a } on the objects with 

S k [{3 a }} = \\ dxdx' £j;(x)0 o (x,x') J^(x'), (2) 

a/3 

where <? (x,x') = G (x,x') - ^ V ® V'G (x,x') is the 
tensor Green function. Next we must constrain the cur- 
rents to be induced sources that depend on shape and 
material of the objects. Formally this is achieved by in- 
tegrating over currents, inserting constraints to ensure 
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that the currents in vacuum simulate the correct induc- 
tion of microscopic polarization P Q and magnetization 
M a (from all multipoles) inside the dielectric objects in 
response to an incident wave. 

Let us consider one object. First, the induced current 
is 3 a = —ikP a + V x M Q , and since P Q = (e a — 1)E, 
M Q = (1 — l//i Q )B, it can be expressed in terms of the 
total fields E, B inside the object as 



-*fc(e Q - 1)E + V x [(1 - l//i a )B] 



(3) 



Second, the total field inside the object must consist 
of the field generated by 3 a and the incident field 
Eo({J a , S a }, x) that has to impinge on the object to in- 



duce 3 a , so that 
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E(x) =E ({J Q ,S a },x)+ifc J dx'£ (x,x')J a (x 

The incident field depends on the current density to 
be induced and on the scattering matrix §" of the ob- 
ject, which connects the incident wave to the scattered 
wave. It is fully specified by the multipole moments 
of 3 a (see below for details). Substituting Eq. Q and 
B = (l/ife)V x E into Eq. (FJ]) yields a self-consistency 
condition that constrains the current J Q . If one writes 
this condition as C Q [J Q ] = for each object, the func- 
tional integration over the currents constrained this way 
for all objects yields the partition function 

Z= f\{V3 a \{ 6(C a {3 a (x)})exp(iS{{3 a }}) . (5) 
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It is instructive to look at two compact objects at a 
distance L, measured between the (arbitrary) origins O a 
inside the objects. In this case the action of Eq. (FJj) is 

Sk[{3 a }] = \J2 J d^ a 3* a ^ a )^-E p ( Xa - L a i) (6) 
+ \Y1 f dx o t dx' a 3a(x o ,)g (x a ,x' a ) Ja(x^) , 

a 

where we have substituted the electric field E a (x ) = 
ik J dx.' a Q (x a ,x.' a ) 3 a (x' a ) and the fields are measured 
now in local coordinates so that x = O a +x a , and L a = L 
(-L) for a — 1(2). The off-diagonal terms in Eq. ^ 
represent the interaction between the currents on the two 
materials. A natural way to decompose the interaction 
between charges is to use the multipole expansion. For 
each body we define magnetic and electric multipoles as 
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- dn a J Q (x Q )V x [n a ji(kr a )Yi* m (i a )} (7) 
(x Q )Vx V x [x a j;(fcr Q )F4(x Q )] , 



for I > 1, |m| < I, where A = ^1(1+1), ji are spher- 
ical Bessel functions and Yj TO spherical harmonics. We 



change variables from currents to multipoles in the func- 
tional integral and, as the final step in our quantization, 
integrate over all multipole fluctuations on the two ob- 
jects weighted by the effective action, 

S k S [{Qf m }} = ~ T U lml> m > Ql'm' ( 8 ) 



Im I'm' 



Qlm ^hnl'm' Ql'm' + X] ("^limVm' Qf'm'} 



with Qf m = (Qujm-)Q%i m )- Let us discuss the terms 
appearing in Eq. Q and sketch its derivation. 

Off-diagonal terms — We need to know the electric 
fields in Eq. ^ exterior to the source that generates 
them. They can be represented in terms of the multipoles 
as E^x,,) = -kT,i m QL^lm(^) where *^(x ? ) are 
outgoing vector solutions of the Helmholtz equation in 
the coordinates of object (3 [12 . We would like to ex- 
press the currents J* in Eq. ^ also in terms of mul- 
tipoles. The difficulty in doing so is that the electric 
field is expressed in terms of outgoing partial waves in 
the coordinates of object /?, while according to Eq. (|7| 



the multipoles involve partial waves $f?(x„) that are 
regular at the origin O a , in the coordinates of object a 
[12] . Going from the outgoing to the regular vector so- 
lutions and changing the coordinate system involves a 
translation and change of basis which can be expressed 
as *^(x a ± Li) = Hvm'Upm'lm^Vm'i^) where the 
universal (shape and material independent) matrices U + 
and U~ represent the interaction between the multipoles. 
For fixed (lm), (I'm 1 ), they are 2x2 matrices (mag- 
netic and electric multipoles), and functions of kL only. 
Their explicit form is known but not provided here to 
save space [13]; they fall off with kL according to clas- 
sical expectations for the EM field. Then the electric 
field becomes ^E^(x Q ± Lz) = Ylim ^L*lm( x «) with 

<t>lm = * Sim ULem'Ql'm" and the integration in Eq. d§| 
leads, using Eq. (Fn), to the off-diagonal terms in Eq. pF. 

Diagonal terms — The self-action, given by the sec- 
ond term of Eq. ^ , is more interesting and more chal- 
lenging. It can be expressed in terms of multipoles if 
we use the constraint for the currents, Eqs. ^ and Q. 
To do so, we first note that in scattering theory one 
usually knows the incident solution and would like to 
find the outgoing scattered solution. They are related 
by the S'-matrix. Here the situation is slightly differ- 
ent. We seek to relate a regular solution E (x Q ) = 
ikJ^im ^o,im*ta(x a ) and the outgoing scattered solu- 
tion, E Q (x Q ) = -kJ2 im Qt m ^° m \x a ), generated by the 
currents in the material — a relation determined by the 
T-matrix, T a = (§ a -I)/2 — schematically iQ a = T Q O 
[HI [15] , We face the inverse problem of determining 

<t>0M 



for known scattering data Qf m , hence 



'JQ,lm 



= ^ Ta \TmVm'Qv 



I'm' 



(9) 
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so that the incident field is given in terms of the S-matrix, 
as indicated in Eq. Q. Next, we express the self-action of 
the currents inside a body (the second term of Eq. pj), 
as S%[J a ] = | J dx a [ED* - BH* - (E D* - B H§)]> 
the change of the field action that results from placing 
the body into the fixed (regular) incident field Eo = Do, 
Ho = Bo, where E, H and D, B are the new total fields 
and fluxes in the presence of the body. Using D = e Q E, 
H = ^q 1 B inside the body and Eq. pj, straightforward 
manipulations lead to the simple self-action 5^[J Q ] = 
— 2lfc / ^ x aJaEo({Ja,§ Q })- If we substitute the regular 
wave expansion for Eo with coefficients of Eq. ([9| and 
integrate by using Eq. ^ , we get Eq. Q . 

The T-matrix can be obtained for dielectric objects of 
arbitrary shape by integrating the standard vector solu- 
tions of the Hclmholtz equation in dielectric media over 
the object's surface |15| and both analytical and numer- 
ical results are available for many shapes [IB] . Hence, 
for the time being, we shall assume that the elements of 
the T-matrix are available. The functional integral over 
multipoles is Gaussian. The resulting partition function 
is an integral over all frequencies of the determinant of 
a matrix M, with inverse T-matrices along the diagonal 
and the matrices U off the diagonal. For each (Im), M 
is a 4 x 4 matrix (2 polarizations for 2 objects). The 
generalization to more than two objects is straightfor- 
ward. The result is formally infinite but the infinity can 
be trivially removed by dividing by Z^, the partition 
function with all objects removed to infinite separations, 
corresponding to setting the off-diagonal terms to zero. 
Dividing by Z^ also cancels the functional Jacobian nec- 
essary to transform from an integral over sources to an 
integral over multipoles. After a Wick rotation, k — > in, 



we finally get the Casimir energy 
tic r °° 



£ = 



2tt 



(iKlogdet^-U^^U+T 1 ) (10) 



in terms of the matrices introduced in Eq. (|8]) . The de- 
pendence of the interaction on distance is completely con- 
tained in V ± , whereas all shape and material dependence 



With 



U- T 2 U+T 1 it 



comes from the T-matrices. 
can be written as £ = f n °° d«TrV°° , ±W which 

ATT JO L^dp— 1 p 

allows for a simple physical interpretation. The matrix 
N scales with distance L as ~ exp(— 2Lk) and describes 
a wave that travels from one object to the other and 
back, involving one scattering at each object. Hence, 
we have obtained a multiple-scattering expansion where 
each elementary two-scattering process, described by N, 
is further decomposed into partial waves. This structure 
allows for a systematic and exact expansion of the in- 
teraction in the inverse distance. At large distance, the 
interaction is determined by the small k scaling of the 



T-matrix, T^ mVm , 



K 



l+l' + l 



This shows that 2p scatter- 
ings become important at order L~ 1 ~~ Gp , and that partial 
waves of order I have to be considered at order L~ 5 ~ 21 . 
Hence, in actual computations, the sum over reflections 
can be cut off at finite p and the matrix N can be trun- 
cated to have dimension 2/(2 + I) x 2/(2 + I) at partial 
wave order I (see below) . We note that Eq. ( 10 1 applies 



also to spatially varying but local e a and fx a , since this 
affects only the T-matrix. Likewise, it can be extended to 
any other boundary conditions or materials by inserting 
the appropriate T-matrix. 

As a specific example, we consider two identical di- 
electric spheres. Due to symmetry, the multipoles are 
decoupled so that the T-matrix is diagonal, 



T, 



I nil). 





Ii+i(nz) +2nzl' l+ i(nz) 


-nl l+ x{nz) 


'l l+h {z)+2zl' l+1 _{z) 




r]K l+ i(z) 


h+ \ ( nz ) + ^ nz I'i+ 1 ( nz ) 


-nl l+ i(nz) 


K l+h (z) + 2zK[ +h {z) 



(11) 



where the sphere radius is R, 



kR, 



V = V e ( iK )/K iK ), and I l+ i, K l+ i 



n = 
are 



y/e(iK)n(iK), 
Bessel func- 



tions. Tj^ lm is obtained from Eq. (11 1 by interchanging 



e and /i. For all partial waves, the leading low frequency 
contribution is determined by the static electric multi- 
pole polarizability, af = [(e - l)/(e + (J + l)//)]i? 2(+1 , 
and the corresponding magnetic polarizability, af 1 — 
[{n-l)/{[i+{l + l)/l)]R 2l+1 . Including the next to lead- 
ing terms, the T-matrix has the structure 



T, 



hnl i 



l)a 



1(21 + iy.\(2i - iy.\ 



and T{% lm is obtained by af 1 



first terms are 7^ 
7& = (4/9)[(m 



7^ 



The 

3 = -[4 + fan + /x - 6)]/[5( M + 2) 2 ]i? 5 , 
l)/(^+2)] 2 i? 6 , and 7^, jf 4 are obtained 
again by the replacement, fj, — > e. Now we can apply our 



general formula in Eq. ( 10 1 to two dielectric spheres with 



r 



center-to-center distance L. For simplicity, we restrict 
to two partial waves (I = 2) and two scatterings (p — 
1), which yields the exact Casimir energy to order L~ 10 . 
Matrix operations are performed with Mathematica, and 
we find the interaction 
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((«l) + («1 ) ) ~ o a i a 



1 



16 L 
315 
'Iff 



(59of - Ha" 



86 7l E 3 



547ft) 



M 



[a? (77 



57JS) 



e «-> m] ■ 



1 

(12) 



where E «-> M indicates terms with exchanged super- 
scripts. The leading term, ~ L~ 7 , has precisely the form 
of the Casimir- Polder force between two atoms [3], in- 



4 



eluding magnetic effects [5J. The higher order terms are 
new, and provide the first systematic result for dielectrics 
with strong curvature. There is no ~ l/L 8 term. 

The limit of perfect metals follows for e — > oo, [i — > 
0. Then higher orders are easily included, yielding an 
asymptotic series 



£ 



hcR 6 



(13) 



where the first 10 coefficients are Co = 143/16, c\ = 
0, c 2 = 7947/160, c 3 = 2065/32, c 4 = 27705347/100800, 
c 5 = -55251/64, c 6 = 1373212550401/144506880, c 7 = 
-7583389/320, c 8 = -2516749144274023/44508119040, 
c 9 = 274953589659739/275251200. This series is ob- 
tained by expanding in powers of N and frequency k, 
and does not converge for any fixed R/L. To obtain the 
energy at all separations, one has to compute Eq. (10) 
without these expansions. This is done by truncating 
the matrix M at a finite multipole order I, and comput- 
ing the determinant and the integral numerically. The 
result is shown in Fig. [I] for perfect metal spheres. Our 
data indicate that the energy converges as e - s ( L / R - 2 ) 1 
to its exact value at I — > oo, with 8 ~ 0(1). Our re- 
sult spans all separations between the Casimir-Poldcr 
limit for L R, and the proximity force approximation 
(PFA) for R/L — > 1/2. At a surface-to-surface distance 
d — AR/'d (R/L — 0.3), PFA overestimates the energy 
by a factor of 10. Including up to I = 32 and extrap- 
olating based on the exponential fit, we can accurately 
determine the Casimir energy down to R/L — 0.49, i.e. 
d = 0.04i?. A similar numerical evaluation can be also 
applied to dielectrics [ID] . 

We have developed a systematic method for comput- 
ing the EM Casimir interaction between compact dielec- 
tric objects of arbitrary shapes. Casimir interactions are 
completely characterized by the S-matrices of the indi- 
vidual bodies. We have computed the force between 
spheres for arbitrary separations, generalizing previous 
results that applied only in singular limits. Our method 
allows for the first time a description of the Casimir inter- 
action from atomic-scale particles (Casimir-Polder limit) 
up to macroscopic objects at short separations (PFA 
limit). For more complicated shapes and multiple ob- 
jects, it would be interesting to probe the dependence 
on the relative orientations of non-spherical objects and 
corrections to pair-wise additivity. Our approach can be 
applied at finite temperatures and extended to the com- 
putation of correlation functions, energy densities, and 
the density of states and may prove also useful to obtain 
thermal (classical) fluctuation forces. 
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FIG. 1: Casimir energy of two metal spheres, divided by the 
PFA estimate £ PFA = -(tt 3 /lU0)hcR/(L-2R) 2 , which holds 
only in the limit R/L — * 1/2. The label I denotes the multi- 
pole order of truncation. The curves / = oo are obtained by 
extrapolation. The Casimir-Polder curve is the leading term 
of Eq. (131. Inset: Convergence at short separations. 
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